Generating optimal states for a homodyne Bell test 
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We present a protocol that produces a conditionally prepared state that can be used for a Bell 
test based on homodyne detection. Based on the results of Munro [PRA 1999], the state is near- 
optimal for Bell-inequality violations based on quadrature-phase homodyne measurements that 
use correlated photon-number states. The scheme utilizes the Gaussian entanglement distillation 
protocol of Eisert et.al. [Annals of Phys. 2004] and uses only beam splitters and photodetection 
to conditionally prepare a non-Gaussian state from a source of two-mode squeezed states with low 
squeezing parameter, permitting a loophole-free test of Bell inequalities. 
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Bell's theorem is regarded by some as one of the most 
profound discoveries of science in the twentieth century. 
Not only does it provide a quantifiable measure of cor- 
relations stronger than any allowed classically, which is 
a key resource in many quantum information processing 
applications, it also addresses fundamental questions in 
the foundations of quantum mechanics. In 1964, Bell 
quantified Bohm's version of the Einstein, Podolsky, and 
Rosen (EPR) gedanken experiment, by introducing an 
inequality that provides a test of local hidden variable 
(LHV) models A violation of Bell's inequality forces 
one to conclude that, contrary to the view held by EPR, 
quantum mechanics can not be both local and real. In or- 
der to experimentally support this conclusion in a strict 
sense, a Bell test that is free from loopholes is required. 
Although it is still quite remarkable that such seemingly 
metaphysical questions can even be put to the test in 
the laboratory, a loophole-free Bell test has yet to be 
achieved. 

For more than three decades, numerous experiments 
have confirmed the predictions of the quantum theory, 
thereby disproving local realistic models as providing a 
correct description of physical reality [2|. However, all ex- 
periments performed to date suffer from at least one of 
the two primary loopholes - the detection loophole and 
the locality loophole. The detection loophole arises due 
to low detector efficiencies that may not permit an ade- 
quate sampling of the ensemble space while the locality 
loophole suggests that component parts of the experi- 
mental apparatus that are not space-like separated could 
influence each other. The majority of Bell tests have used 
optical systems to measure correlations, some achieving 
space-like separations but still subjected to low efficiency 
photodetectors (see, e.g., Ref. §)• Correlations in the 
properties of entangled ions were shown to violate a Bell 
inequality using high efficiency detectors eliminating the 
detection loophole; however, the ions were not space- 



like separated jj]. A major challenge that has yet to 
be achieved is to experimentally realize a single Bell test 
that closes these loopholes. 

The ease with which optical setups address the locality 
loophole coupled with the currently achievable high effi- 
ciencies (> 0.95) of homodyne detectors make Bell tests 
using quadrature-phase measurements good candidates 
for a loophole-free experiment. Furthermore, continuous 
quadrature amplitudes are the optical analog of position 
and momentum and more closely resemble the original 
state considered by EPR. Unlike photon counting exper- 
iments which deal with the microscopic resolution of a 
small number of photons, by mixing the signal with a 
strong field, homodyne measurements allow one to de- 
tect a macroscopic current 

In this article, we propose a test of Bell inequalities us- 
ing homodyne detection on a conditional non-Gaussian 
"source" state, prepared using only passive optics and 
photon detection. Events are pre-selected - using event- 
ready detection one knows with certainty that the de- 
sired source state has been produced - requiring no post- 
processing. Photon detectors are only used in the pre- 
selection process and only affect the probability of suc- 
cessfully creating the source state whereas the actual 
correlation measurements are performed using high ef- 
ficiency homodyne detectors. The source is a correlated 
photon-number state that is near-optimal for Bell tests 
using homodyne detection, opening the possibility of a 
conclusive, loophole-free test. 

We consider a two-mode quantum state of light that 
can be written as 



I*) = ^CnK™)' 



(1) 
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which is correlated in photon number \n,n) = \n) A®\n) b 
for modes A and B. For example, the two-mode squeezed 
state \ip\) has coefficients given by c„ = A n \/1 — A 2 , 
where A = tanh(s) is determined by the squeezing pa- 
rameter s @ . Such states are experimentally easy to gen- 
erate; however, because they possess a Gaussian Wigner 
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distribution in phase space, they are unsuitable for tests 
of Bell inequalities using quadrature-phase measurements 
as it is a requirement that the Wigner function pos- 
sesses negative regions pj. Alternative, theoretically pre- 
dicted two-mode quantum superposition states called cir- 
cle states, also generated from vacuum fields through 
nondegenerate parametric oscillation, having coefficients 
given by c„ = r 2n / 'nl-J 'Io(2r 2 ) , do exhibit a violation 
for quadrature-phase measurements with a maximum vi- 
olation occurring for r = 1.12 Unfortunately, un- 
like the two-mode squeezed states, circle states are dif- 
ficult to realize experimentally. A recently proposed so- 
lution towards an experimentally realizable state of light 
that is suitable for a homodyne Bell test is the photon- 
subtracted two-mode squeezed state || Q , having coeffi- 
cients c„ = ^/(l - A 2 ) 3 /(l + A 2 )(n+ 1)A", which utilizes 
non-Gaussian operations on a Gaussian state. In this 
scheme, a photon is detected from each mode of a two- 
mode squeezed state and only the resulting conditional 
state is used for correlation measurements in the Bell 
test. While the two-mode squeezed state has a positive- 
everywhere Wigner function, the conditional state after 
photon subtraction does not. 

To date, all proposed states for a Bell test using 
quadrature-phase measurements are not optimal states, 
meaning that they do not produce the maximum possi- 
ble violation of Bell inequalities. The scheme presented 
here produces a two-mode photon-entangled state that 
is near-optimal, using only beam splitters and photon 
detection. The beam splitter may be described by the 
unitary operator [Io| 

tt rria* a — R* a Rba*rr>—b*b /o\ 

U a b — J- e el, (2) 

which describes the mixing of two modes a and b at 
a beam splitter with transmissivity T and reflectivity 
R. On-off photon detection is described by the positive 
operator- valued measure (POVM) of each detector, given 
by 

n = |o)(o|, n 1 = 7-|o)(o|. (3) 

The on-off detectors distinguish between vacuum and the 
presence of any number of photons. The procedure is 
event-ready, a term introduced by Bell, in the sense that 
one has a classical signal indicating whether a measurable 
system has been produced. The states demonstrating a 
violation of local realism presented here do not rely on 
the production of exotic states of light; in fact, only a 
parametric source generating states with a low squeezing 
parameter is required, making the procedure experimen- 
tally feasible with current technology. As depicted by 
Fig. 1, there are three parties involved: Alice, Bob, and 
Sophie. Sophie prepares the source states that are sent to 
Alice and Bob, who perform correlation measurements. 
We first describe the procedure Sophie uses to generate 
the source states, which is shown by the diagram in Fig. 
2, and then discuss the measurements performed by Alice 
and Bob. 
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FIG. 1: A schematic diagram of the conditional homodyne 
detection for a Bell test. Sophie records successful prepara- 
tion of the source state, denoted by \i>^)ps, while Alice and 
Bob perform homodyne measurements at their space-like sep- 
arated locations. The three parties agree on which events to 
discard based on their measurement records. 

In the first step, two-mode squeezed states are mixed 
pairwise at unbalanced beam splitters followed by the 
non-Gaussian operation associated with the POVM ele- 
ment III. Specifically, a non-Gaussian state is generated 
by 

(n 1)C ® n M )(t/ ac ® UuWx) IV>a>, (4) 

where \ip\) denotes the two-mode squeezed state with 
c„ = A n vl — A 2 . For sufficiently small A, the operator 
III describing the presence of photons at the detector ap- 
proaches the rank-one projection onto the single photon 
number subspace |1)(1|, which is still a non-Gaussian op- 
eration. Under this condition, (un-normalized) states of 
the form 

|^ 0) ) = |0,0)+e|l,l) (5) 

can be produced. That is, even though the output state 
of (0J will in general be a mixed state, when A £ [0, 1) is 
very small, the resulting states can be made arbitrarily 
close in trace-norm to an entangled state with state vec- 
tor given by Eq. JSJ, provided the appropriate choice of 
beam splitter transmittivity |T(A)| = |£- v^ 2 + 8A 2 |/4A 
is used llj. It should be emphasized that the state ), 
having a Bell-state form, can be generated for arbitrary 

It is interesting to note that the state given by Eq. (J5J 
does not violate a Bell inequality for quadrature-phase 
measurements for any £, even when it has the form of 
a maximally entangled Bell state, as was shown in Ref. 
|13|. in which a numerical study of the optimal coeffi- 
cients for Eq. (JIJ was performed. For certain values 
of £, Eq. (0) describes a state that possesses a Wigner 
distribution that has negative regions, showing that neg- 
ativity of the Wigner function is a necessary but not suf- 
ficient condition for a violation of Bell inequalities using 
quadrature-phase measurements. 

The second step is to combine two copies of the state 
given by Eq. ||SJ pairwise and locally at 50:50 beam split- 
ters described by the unitary operator of Eq. (0). De- 
tectors that distinguish only between the absence and 
presence of photons are placed at the output port of 
each beam splitter and when no photons are detected, 
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the state is retained. The resulting un-normalized state 
is 

oo 

= (o,o|c/ ac ® tf M |V w >|V w > = E c « +1) m, 

n=0 

(6) 

where the coefficients are given by [14| 

C , ( i i+1) =2-"E(;)c«cll. (7) 

It is optimal to iterate this procedure three times so that 
Sophie prepares the state \4>^ 3 ')- Each iteration leads 
to a Gaussification of the initial state [Til IT^ . which 
builds up correlated photon number pairs in the sum of 
Eq. QJ. Further iterations would Gaussify the state too 
much and destroy the nonlocal features for phase space 
measurements . 

The final step is to reduce the vacuum contribution by 
subtracting a photon from each mode of the state IV''- 3 ''), 
obtaining a state proportional to ab\ip^). This is done 
by mixing each mode with vacuum at two beam splitters 
with low reflectivity. A very low reflectivity results in 
single photon counts at each detector with a high prob- 
ability when a detection event has occurred. Thus, the 
unitary operation describing the action of the beam split- 
ter is expanded to second order in the reflectivity and the 
state is conditioned on the result N = 1 at each detector. 
The final photon-subtracted state, given by 

oo 

|^ 3 ))p5=^E(n+l) C i 3 j 1 |n,n), (8) 

where J\f is a normalization factor, is a near-optimal state 
for homodyne detection. Figure 3 compares the previ- 
ously proposed states - the circle state and the photon- 
subtracted two-mode squeezed state - with the near- 
optimal state |?/>( 3 )) ps, as well as the numerically op- 
timized state in Ref. |l3j . The conditioning procedure 
alters the photon number distribution of the input state 
and behaves similarly to entanglement distillation. 

Although the procedure used to create the correlated 
photon source is probabilistic, with the success proba- 
bility determined by the amount of two-mode squeez- 
ing and the transmittivity of the unbalanced beam split- 
ters, it is event-ready - Sophie has a record of when the 
source state was successfully prepared. Low efficiency 
photon detectors used in the state preparation only af- 
fect the success probability and do not constitute a de- 
tection loophole. Each mode of the source state \ip^ 3 ')ps 
is distributed to a separate location where correlation 
measurements using high efficiency homodyne detectors 
are performed by the two distant (space-like separated) 
parties, Alice and Bob. Alice and Bob each mix their 
light modes with independent local oscillators (LO) and 
randomly measure the relative phase between the beam 
and the LO, taking into account the timing constraint 
that ensures fair sampling. Alice measures the rotated 




FIG. 2: A tree diagram of the conditional preparation of 
the state denoted by from a finite supply of two-mode 

squeezed states |i/>^). Each circle represents the two-mode 
state as given on the right. The black (white) circles cor- 
respond to the Gaussian (non-Gaussian) operations used to 
produce the states at each node. 

quadrature xf = x A cos 9 + p A sin 9 and Bob measures 
the rotated quadrature x B = x B cos (j> + p B sin cf>. Corre- 
lations are considered for two choices of relative phase: 
9\ or 02 for Alice and </>i or <p2 for Bob. Finally, Alice, 
Bob and Sophie compare their experimental results to 
determine when the source state was successfully gener- 
ated and which correlation measurements to use for the 
Bell inequalities. 

Two types of Bell inequalities will be examined - 
the Clauser-Horne-Shimony-Holt (CHSH) and Clauser- 
Horne (CH) inequalities [Ta|. To apply these inequalities, 
which arc for dichotomous variables, the measurement 
outcomes for Alice and Bob are discretized by assigning 
the value +1 if x > and -1 if x < 0. Let P^f (#,</>) 
denote the joint probability that Alice and Bob realize 
the value +1 upon measuring 9 and 4>, respectively and 
P A {9) denote the probability that Alice realizes the value 
+1 regardless of Bob's outcome, with similar notation 
for the remaining possible outcomes. From LHV theo- 
ries, the following joint probability distribution can be 
derived: 

P AB (9,<f>) = J p(X)pf (9, X) P f ((f), X)d\ (9) 

with i,j = ±, by postulating the existence of hidden vari- 
ables A and independence of outcomes for Alice and Bob. 
Quantum mechanically, the joint probability distribution 
is given by the Born rule P(x A ,x B ) = \(x A , x B \i/j^ 3 ')ps\ 2 - 
The probability for Alice and Bob to both obtain the 
value +1 is P AB (9, 0) = / °° J °° P(xf, x B )dx A dx B . The 
joint distribution is symmetric and a function of only the 
sum of the angles x — @ + 4> permitting the identifica- 
tion P AB (9^) = P AB ( X ) = P AB (~ X ) and P AB ( X ) = 
P AB (x)- The marginal distributions Pf(0) = Pf (</>) = 
1/2 are independent of the angle. Given the probability 
distributions, the predictions of quantum theory can be 
tested with those of LHV theory. 

First, we consider the Bell inequality of the CHSH 
type, which arises from linear combination of correlation 
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functions having the form 

B = E{9 u 4n) + E{6x,<t> 2 ) + E(e 2 ,4>i) - E(6 2 , <fe), (10) 

where E(0i, <f)j) is the correlation function for Alice mea- 
suring 6i and Bob measuring <pj. These correlations are 
in turn determined by 



E(0,<t>) = P 



-P 



AB i 



V ; 



)-P±*(6,, 



(11) 

obtained through the many measurements that infer the 
distributions P^ B (9,(f>). With the aid of the symmetry 
and angle factorization properties, the CHSH inequality 
takes the simple form B = 3E(x) — E(3x) with LHV 
models demanding that \B\ < 2. The strongest violation 
of the inequality is obtained for the value \ — 7r /4, thus, 
a good choice of relative phases for Alice and Bob's mea- 
surements is Q\ = 0, 6 2 = 7r/2, (pi = — 7r/4, and <fi 2 = 7r/4. 
Using homodyne detection with optimal correlated pho- 
ton number states, the maximum achievable violation is 
2.076 whereas using the source states presented here, a 
Bell inequality violation of B = 2.071 is achievable. 

Let us also consider the Clauser-Horne (strong) Bell 
inequality formed by the linear combination 



l) 



h)+p$n 



Pf(0 2 )+Pf(^l) 

(12) 

denoted by S, for which local realism imposes the bound 
|5| < 1. Again, using the properties of the probability 
distributions, the simplification S = (x) — P++ (3x) 

is possible. With the following choice of the phases: d\ — 
0, 9 2 = 7r/2, (j>i — — 7r/4, and 4> 2 = 7r/4, a violation 
of S = 1.018 is attainable given the states in Eq. (J5J 
with parameter value £ = l/y/2, which is quite close to 
the maximum value of 1.019 achieved by the numerical, 
optimal states in Ref. [LJ • 

We have shown how it is possible to prepare a near- 
optimal state for a Bell test that uses quadrature-phase 
homodyne measurements. Only very low squeezed states, 
passive optical elements and photon detectors are re- 
quired, making the procedure experimentally feasible 
at present. An initial state with a positive-every where 
Wigner function was succeeded by both non-Gaussian 



and Gaussian operations to prepare a state that exhibits 
a strong violation of both the CHSH and CH Bell in- 
equalities. Efforts are currently being made towards an 
experimental realization of entanglement distillation for 
Gaussian states. The procedure presented here offers the 
opportunity for another possible experiment, as it uti- 
lizes a subset of an entanglement distillation procedure. 
Of course, any observed violation of a Bell inequality 
is sensitive to inefficiencies in the experiment that tend 
to deplete correlations. A full analysis involving dark 
counts and detection inefficiencies as addressed in Ref. 

is necessary. Near-optimal states for homodyne de- 
tection may allow a larger window for experimental im- 
perfections and offer the opportunity for a conclusive, 
loophole-free Bell test. 
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FIG. 3: The coefficients in Eq. Q are plotted as a function 
of photon number n. The curves represent photon-subtracted 
state two-mode squeezed state (small dashes) for A = 0.6, the 
circle states (large dashes) for r — 1.12, and the two-mode 
squeezed state (dashed-dotted) for A = 0.6. The optimal 
states are shown by squares and the source states |^ 3 ')ps 
are shown by grey triangles for £ = 0.71. 
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